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Abstract 

The covariant quark model represents an appropriate theoretical framework to describe the recent results 
onB s — > J/ip + n and B s — > J/ip + n' decays from the Belle collaboration. In this article we present the main 
features of the covariant quark model together with details on some of its aspects and methods, which we 
consider to be important. Further we apply the model specifically to the studied decay processes and give 
numerical results on decay widths as they follow from the model. We conclude that the model, with most of 
its parameters previously fixed from different processes, is able to incorporate the new measurements from 
the Belle experiment. 
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I. INTRODUCTION 



The low energy region of quantum chromodynamics (QCD) is an interesting and challenging 
area of investigation nowadays. On one hand many accurate data exists in this domain, especially 
results from hadron spectroscopy and hadron decays, which might give a handle to understanding 
the quark dynamics inside hadrons. Presently the precision and the amount of these data con- 
tinuously increases thanks to ongoing high-energy experiments and heavy-quark factories. On the 
other hand theoretical predictions encounter difficulties: since with decreasing energy the coupling 
constant increases, the perturbative approach (pQCD) looses its applicability and cannot provide 
valid results. Also lattice calculations are generally seen as not yet enough developed and precise to 
be considered as a well established low-energy solution of the QCD which would satisfactory explain 
its numerous features. Undoubtedly however, some theoretical description of hadron dynamics is 
needed at least for a correct description of physic and background measured in particle detectors. 
One thus has to rely on a model-dependent approach. In this text we use the covariant quark model 
to study the decays B s — > J/ip + rj and B s — > J/ip + r]', processes that have been recently measured 
at Belle [l|. 

The covariant quark model jgj] is an ambitious model with many appealing features. Its 
Lagrangian-based formulation leads to a full Lorentz invariance and, in the limit of large num- 
ber of hadrons, it has only one free parameter per hadron. The model has wide applications 
spectra and was already successfully used to describe different types of heavy hadron decays and 
other hadron-related observables (magnetic moments, cross sections, tetraquark properties, etc..) 
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II. COVARIANT QUARK MODEL 

The interaction Lagrangian (density) of the model 

£ int = g H H (x) • J H (x) + H.c. (1) 
is constructed from the hadron field H (x) and the quark current. The latter is in case of mesons 



27l | written as 

Jh(x) = j dxi j dx 2 F M (x;x 1 ,x 2 ) • q^ (x x ) r M q| (x 2 ) (2) 

and makes appear that, within the model, the interaction is mediated only by the quarks and the 
gluons are absent. The form of the vertex function Fm is chosen such as to reflect the intuitive 



expectations about relative quark-hadron positions 

F M (x;xi,x 2 ) = 5 (4) (x - wixi - w 2 x 2 ) $ M (xi - x 2 ) 2 , (3) 
where we require wi + w 2 = 1 . We actually adopt the most natural choice 

111; . , 

Wi = -!— , i = 1,2 4 

mi + m 2 

where the barycenter of the hadron is identified with the barycenter of the quark system. The 



(xi - x 2 ) 2 



is assumed to have a Gaussian form which is in the momentum 



interaction strength <I>m 
representation written as 

<I M (-P 2 )=exp(i0. (5) 

Here Am is a hadron-related size parameter which is regarded as a free parameter of the model. 
Additional free parameters are the quark masses 

m Ujd = 0.235 GeV, m s = 0.424 GeV, 

m c = 2.16 GeV, m b = 5.09 GeV (6) 

and a universal cut-off parameter A cu t- ff , which is commented on in more details later in this text. 
The numerical values of these parameters as well as of size parameters for several hadrons were fixed 
by fitting the model to measurement data of simple quantities, such as leptonic decay constants and 
electromagnetic decay widths fl^ . The value of size parameter for B s meson is Ab s = 1.95 GeV. 
The model thus has in total Nh + 5 parameters: for each of Nh hadrons one A parameter, four quark 
masses and one universal cut-off. The coupling constants gM can be related to these parameters 
using the so-called compositeness condition, which is discussed in the following section. 



III. COMPOSITENESS CONDITION, CALCULATION METHODS AND INFRARED 
CONFINEMENT 

In this section we give more details about some issues related to the model. As announced 
previously, one of them is the compositeness condition. 

The quark fields as well as the hadron field enter the interaction Lagrangian of the model ([T|) 
as elementary although, in nature, the hadrons are compound of quarks. An effort was made in 
theoretical physics to find an appropriate description of composite particles. The authors of the 
references [la . \l\ argue, that the hadron fields renormalization constant can be interpreted as 
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the matrix element between the physical state and the corresponding bare state. The case Zh = 
thus corresponds to a state not containing the bare state and can be therefore properly interpreted 
as a bound state. This idea was used and introduced into the covariant quark model 28|. The 



renormalization constant is expressed through the derivative of the meson mass operator and takes 
the form 

Zm = 1 - §n M W) = 0. (7) 

The coupling constants are in this way eliminated as free parameters. This not only gives the model 
more predictive power but also helps to stabilize model predictions over wide spectra of hadron data. 
The first step in calculation of physical observables is thus the determination of couplings gM of 
participating hadrons. 

In computation of Feynman diagrams we use, in the momentum space, the Schwinger represen- 
tation of the quark propagator 

oo 

= (m + k) f dae- a ( m2 ~ k2 ) (8) 
o 

The general form of a resulting the Feynman diagrams is 

n( Pl ,...,p m ) = j d n aj [d^j^Xexpj-^aiJmf-^i + Pi) 2 ]!, (9) 

where Kj represents a linear combination of loop momenta, P; stands for a linear combination of 
external momenta and <3? refers to the numerator product of propagators and vertex functions. The 
evaluation of these expressions can be much simplified if done in a smart way. One can use two 
operator identities, of which the first one 



j d 4 kP(k) e 2kr =y d 4 kP 



1 9 \ okr 

2<9r ' 



2dr 

is suited for en elegant loop momenta integration. The second one 



P(^) /d 4 ke 2 ^ (JO. 



oo 



d n « e-4 == I d«a e~4p (11 - I ] 1, ( J ., , 



2 di J J V 2 dv a 

o o 

where r = r (a;) and a = a (Am, simplifies the computation following the trace evaluation: the 
polynomial in the derivative operator which results from the trace can be applied to an identity, 
instead being applied to a more complicated exponential function. 



The last point which remains to be discussed is the cut-off we apply in the integration over the 
Schwinger parameters. This integration is multidimensional with the limits from to +00. In order 
to arrive to a single cut-off parameter we firstly transform the integral over an infinite space into an 
integral over a simplex convoluted with only one-dimensional improper integral. For that purpose 
we use the ^-function form of the identity 

1 = J dtS {t - Yj a ^j ( 12 ) 

from which follows 

11 = j dtt 11 - 1 j d^aslt-J^aA xW (tax,..., tax), (13) 
\ 1=1 J 

where W represents the integrand of Schwinger parameters. The cut-off A is then introduced in a 

natural way 

00 1/A 2 


Such a cut-off makes the integral to be an analytic function without any singularities. In this way 
all potential thresholds in the quark loop diagrams are removed together with corresponding branch 
points [2|]. Within covariant quark model the cut-off parameter is universal for all processes and its 
value, as obtained from a fit to data, is 

Acut-ofr = 0.181 GeV. (15) 
The integrals are computed numerically. 



IV. DECAYS B s -> J/ip + T) AND B s -> J/ip + 7/ 

The diagram of the B s — V J/ip + r/W decay is shown in Fig. [Taj The 77 and r( particles are 
mixture of the light and the s-quark component. In the approximation m u = = m q one can 
write their quark content as 

77 : -= sin 5 (uu + dd) — cos 5 (ss) 

= — sin 5 (qq) — cos 5 (ss) , (16) 



cos 5 (uu + dd) — sin 6 (ss) 
= cos 5 (qq) — sin 5 (ss) , 
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(17) 



with 



and 



s s 

(a) (b) 
Figure 1: a) Diagram of B s — V J ftp + i]^ decay, b) Factorization of the diagram. 



qq = ±= (ml + dd) (18) 



6 = e P - e h 

9\ = arctan —= m 35.26°, 

6p& -13.34°. (19) 

The value of the pseudoscalar angle 6p is deduced from the Ref. {is| . where however a different 
convention is used for the mixing angle ipp. It can be related to ours by 5 = ipp — tt/2. 

We describe the decay only through the dominant s-quark contribution to the rf'' meson since 
the light quark one results from higher-order diagrams. The quark-hadron interaction is described 
by the Lagrangians from the Section [XT] (see Appendix [AT]) and we use an effective theory with four- 
quark interaction to describe the production of the J ftp particle. The Lagrangian of this interaction 
is written as 

Aff^VcbV^CiQi, (20) 



where Gf is the Fermi Coup 



)ling Constant, V xy refers to the elements of the CKM matrix, Ci are 
the Wilson coefficients (Ref. 16]) and Qi are operators listed in Appendix I A 21 To each component 
of both r\ and r( particles we associate a size parameter. So, in general, the description requires 
four free parameters A^ q , A^ s , A^? and A^. We consider the mixing angle 5 as fixed although, 
in principle, it also might be varied to determine its most suited value within the covariant quark 
model. As a model-independent parameter it can be compared to other models and to data. 

An important simplification in calculations comes from the evaluation of the Qi-operator matrix 
elements. In can be shown that the expression factorizes into a part that corresponds to the 
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4.67 x 10~ 4 
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4.04 x 10~ 4 
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Table I: Decay widths and branching fractions for selected processes: model predictions and data. 



B s — > 77 W transition form factor and a part that is proportional to the decay constant of the J ftp 
particle (F igllbl and Appendix [A3]). One can so refer to existing formulas for form factors and decay 
constants [12| and re-use previous computer programs for a numerical evaluation with different 
settings. 

In addition to the Belle results, we need further data in order to over-constrain the model. We 
have chosen the following decays 



V 

rf- 
P° 

CO 

ri- 



ll 
11 

• VI 

• VI 



tp — >■ 777 

ip -4 7/7 

B d + 

B d -4 J/V> + rf 



(21) 



These processes have already been previously described by the covariant quark model [2j] and can 
be straightforwardly included into a global fit. 
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V. RESULTS AND DISCUSSION 



The optimal model parameters (in GeV) 

Af = 0.881, = 1.973 

AJ? = 0.257, A^ = 2.797 (22) 

were obtained from a \ 2 fit to the data. The corresponding description of the data by the covariant 
quark model is shown in Table HJ 

The interpretation of the results might by done in several steps. Firstly, one observes that all 
model predictions are of the order of the experimental numbers. Majority of them have the relative 
error smaller then 30%, when compared to the data. The most important relative error is 73% 
in case of T^r^^-j, still significantly smaller then factor two. The latter case actually suggests one 
might want to consider the possibility of a gluonium content of the rf'> meson, as discussed in 
the Ref. [15l | (and other references therein). On the other hand one must admit that, when the 
experimental errors are taken into account, the model prediction are usually quite outside the error 
intervals. Here one can argue in two ways. Firstly, the model is only an approximation to the 
first-principle theory. One thus, from beginning, does not expect the model to be fully accurate. 
Consequently, the goodness of the model when interpreted through the data error intervals depends 
on the data precision measurement. Every approximate model becomes "very bad" when the data 
becomes very precise. So the point of view based on the error intervals might not be the most 
suited one. 

An optional criterion might be a comparison to the experimental needs. Experiments usually 
need a model that allows for an appropriate correction of detector effects. In fact they usually need 
more models so as to be able to establish a systematic error related to the data correction. We think 
that within this logic, the results of the covariant quark model make the model fully acceptable 
and legitimate is development and eventual use. 

Yet, an additional approach to rate the model is to compare it to other existing models. We have 
chosen some Q-Q of not very numerous works, that give explicit numbers for a set of observables 
overlapping with those chosen by us. When comparing processes in common (Table |T|), none of this 
models describes the data better than ours (if total the x 2 is calculated). This fact confirms, that 
the covariant quark model is a very competitive one among the available models. 

Finally, one can still await further confirmation and more precision in the experimental data, 
especially in the case of recent results obtained by a single collaboration which have not yet been 
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independently cross-checked. The final picture concerning the data description by the model might 
still change. 
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Appendix A: Expressions and Formulas 

We make use of this Appendix to display longer formulas referred in the text. 

1. Lagrangian of the model 

The Lagrangian (density) is written as 

£ = £b s + A/ + £j/t/> + £ e fr 



with C e g being given in the text. 



■^0 



(*)// 



£ Bs (x) = gB s B s (x) // dx!dx 2 <5 (x - w b x x - w s x 2 ) 4> Bs (xi - x 2 ) b (xi) i7 5 s (x 2 ) + h. c 



£r) (x) = g v V ( x ) jj dxidx 2 5 ^x - ^xi - ^x 2 ^ (f) v (x x - x 2 ) 



> < --j= sin (5) [u (xi) i7 5 u (x 2 ) + d (xi) i7 5 d (x 2 )] - cos (5) [s (xi) i7 5 s (x 2 )] 



C lf (x) - g,,<// (x - // dxidx 2 <5 ^x - ^xi - ^x 2 ^ (f) v , (xi - x 2 ) 

> <, cos (5) [u (xi) i7 5 u (x 2 ) + d (xi) i7 5 d (x 2 )] - sin (S) [s (xi) i7 5 s (x 2 ) 



'j/V = SV-^m ( x ) dxidx 2< 5 fx - ^xi - ^x 2 



(xi -x 2 ) 2 c(xi)7^c(x 2 ) 
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2. Four-quark vertex operators 

The four-quark operators read as follows 



with 



Qi 


— ( c aib a2 


)v-A (Sa 2 c ai)v-A 


Q 4 = 


{s ai b a2 )y 


_A (Ca 2 C ai )y_ J 4 


Q 2 


= (Caib ai 


)v-A (Sa 2 c a 2 )v-A 


Q 5 = 


(Saib ai )v- 


_A (Ca 2 C a2 )v+A 


Q 3 


= (s ai b ai 


)v-A (Ca 2 C a2 )v-A 


Qe = 


(Sai b a2 )y. 


-A (Ca 2 C ai )y +J 4 


(^)v-a 


= rpO^ip, 


0^ = 7^ (1 - 7 5 ) 


(^)v- 


f A = 


^ 0^ = 7 M (l + 7 5 ) 



3. Some elements on factorization 

For the r\ particle and the operator Qi, the time-product matrix element is written as 

(t { [b (x?) i 7 5 s (xf )] Cl Q! (x-) [g (4) i7 5 s (4)] [e (4) w {4)}}) ■ 

This, after evaluating the contractions and color indices, leads to 

(T{. . . Qi . . .}) = -CiN^,„ x Tr [S b (x w - xf ) 7 5 S S (xf - x?) 7 5 S S (4 - x w ) O"] 



xTr 



S c (x w -xn 7 ly S c (xf -x w )0^ 



where S is a propagator, index w refers to the position of the four-quark interaction and e^ u is 
the J/ip polarization vector. The structure of the expression makes visible the factorization into 



a "B s — rf and u J/tp" part. When actually comparing the two parts to the expressions from 
one recognizes the form factor and the decay constant. Situation is analogical for rf and other 
operators. 



Further, very simple relations exist between the matrix elements for different operators. One 



has 



<T{...Q 2 ,4, 6 ...}) C 



N c 



(T{...Q 1>3) 5...}} . 
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